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Abstract 

New methods for D-decomposition analysis are presented. They are 
based on topology of real algebraic varieties and computational real al- 
gebraic geometry. The estimate of number of root invariant regions for 
polynomial parametric families of polynomial and matrices is given. For 
the case of two parametric family more sharp estimate is proven. The- 
oretic results are supported by various numerical simulations that show 
higher precision of presented methods with respect to traditional ones. 
The presented methods are inherently global and could be applied for 
studying D-decomposition for the space of parameters as a whole instead 
of some prescribed regions. For symbolic computations the Maple v. 14 
software and its package RegularChains are used. 

1 Introduction 

In this paper linear differential equation with constant coefficients 

x(t) = Ax(t), 

where A is a complex(real) nxn matrix, will be called linear continuous system. 
Similarly, linear discrete system is a finitc-diffcrcncc equation 

x(t + 1) = Ax(t), 

where A is a complex(real) nxn matrix. 

The problem of finding stability domains of a linear control system depending 
on a vector parameter is a classical one. It can be formulated as such: Let 
us consider a polynomial P(s, k) with complex(for example, real) coefficients, 
polynomially depending on real parameters k — (k\, . . . ,ki) £ K. Complex 
parameters could be considered as a + jb, where a, b are real. 

This polynomial is usually considered as a characteristic polynomial of a lin- 
ear system. Vector k understood as indeterminacy or parameters of a regulator. 

We need to find such regions U C K that for all k £ U P(s, k) will be 
continuous-time(resp. discrete-time) stable. 

Polynomial is called continuous-time stable(Hurwitz) if all of its roots have 
negative real part. Polynomial is called discrete-time stable(Schur) if all of it's 
roots is have modulus lesser than 1. 
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Continuous-time(res. discrete-time) system is stable if and only if its char- 
acteristic polynomial is a Hurwitz(resp. Schur) one. 

It is obvious that the problems of finding regions of hurwitzness or schurness 
of a polynomial is equivalent. Clearly, polynomial P(z) is Schur if and only if 
(s — l) dc s- p 'p(*±I) is Hurwitz. Therefore we can consider only continuous-time 
systems, and we will do that in any case, unless otherwise specified. 

There are two cases which are of the greater interest for control theoretist, 
namely the case of characteristic polynomial of SISO system with affine depen- 
dency on parameters and the case of family of matrices affinely depending on 
parameters, for example, in the case of construction of a linear output feedback. 
In the latter case polynomial P(s, k), will be det(si — (A + ki^-i))- 

One of a methods of studying stability regions is a D-decomposition ap- 
proach based on study of decomposition of parameter space into connected 
regions with equal numbers of stable and unstable roots(root-invariant regions 
or regions of ^-decomposition) . It is clear that the border of these regions is 
a semialgebraic hypersurface in K parametrized by equalities 5ft(P(jw, k)) — 
0, $$(P(ju!, k)) = and, possibly, a hypersurface with an equation \adegp(k)\ = 
0, where a^cgP - is a leading term of P. 5ft and 9 here denote real and imaginary 
part of an expression. 

The main idea of this approach goes back to I. Vishnegradsky[30 , but first 
clear formulation and detailed study of a method have been done by Yu.I. 
Neimark [2T], [22]. The contemporary state-of-the-art could be obtained from 
paper [T3] and review [T3]. 

The main idea of our paper is an application of various methods of compu- 
tational algebra and real algebraic geometry to the study of D-decomposition- 
both computational(finding sample points from all root invariant regions) and 
theoretical - upper bounds to the number of regions. 

The famous Tarski-Seidenberg theorem ([27] [25]) tells us that the problem of 
description of all root-invariant regions, as well as stable ones is algorithmically 
solvable. Nevertheless, effective algorithms of study of real polynomial systems 
of equations have been included in general purpose computer algebra systems 
only recently. For example, Maple package RegularChains was included in 
Maple only in 2005 and it is still rapidly developing. The history of those 
algorithms itself goes back to discovery of cylindric algebraic decomposition 
by G.E. Collins in 1975 [7] and partial cylindric algebraic decomposition by 
G.E.Collins and H.Hong in 1991 8 . Monograph [3] is a standart reference book 
on the subject. 

Moreover, applications of quantifier elimination algorithms to control theory 
goes back to the seminal paper by B.D.O. Anderson, N.K. Bose, E.I. Jury pub- 
lished in 1975 [2J, where an application of quantifier elimination to the linear 
output feedback stabilisation problem have been studied. They used Routh- 
Hurwitz theory, multivariate resultants and classic Tarski-Seidenberg quantifier 
elimination technique to check existence or non-existence of stable output feed- 
back. However, there was no effective computational techniques, software and 
sufficiently powerful machines at the time, so the interest to this topic have 
risen only in the midde of 1990's[TT] [23] [IB]. The application of real algebraic 
geometry to the study of the stable points set and the Nyquist map is stated 
as an unsolved problem in monograph by E.A. Jonckheere [19] The approach 
used in this papers is quite similar to the Anderson-one, but authors have used 
not an original Tarski-Seidenberg algorithms but cylindrical algebraic decom- 
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position implemented at that time in the QEPCAD software by H.Hong. His 
collaborative work with R. Liska and S.Steinberg [16 gave rise to applications 
of partial cylindric algebraic decomposition algorithms to stability problems for 
PDE. 

Later, H.Anai and S.Hara et al. in series of papers have introduced the 
new methods of robust control synthesis with quantifier elimination based on 
so-called Sturm-Habicht sequence, Symbolic-Numeric CAD etc (see PQ,[H] and 
references therein), they have done computations with these methods in com- 
puter algebra system Risa/Asir. They have also designed a Maple package 
SyNRAC, where some of their algorithms was implemented. 

Some similar considerations in the context of Nyquist criterium could be 
found in papers by N.P. Ke ([50] and references therein). The paper is 
devoted to the applications of computational real algebraic geometry to an 
asymptotic stability analysis of nonlinear systems. A.D. Bruno in his preprint 
[B] introduces use of resultants for study of D-decomposition border. 

This paper is organised as follows. In section 2 we obtain some upper 
bounds to the number of root invariant regions in the parameter space in 
two-dimensional and general cases. Using both Groebner bases and classical 
resultant-based elimination technique we obtain a variety containing projection 
of intersection of hypersurfaces $l(P(juj, k)) = 0, 3(P(ju;, k)) = onto K. Then 
we estimate the number of regions of its complement using classical Harnack and 
Bezout theorems in two-dimensional case and Warren inequality in the general 
case. These results are far-going generalisation of upper bounds from [T2l 115] . 

In section 3 we give the description of Maple-based algorithm of study- 
ing D-decomposition and concentrate on some particular examples from |13) . 
which illustrate main features of our approach. We use a combination of Groeb- 
ner basis-based or resultant-based elimination and partial cylindric algebraic 
decomposition. Specifically, for the continuous-time system first appeared in 
|28j . Example 2 and studied in [T5J, Example 15, we have found one new root- 
invariant region. This example tell us about one very important feature of our 
methodology - its inherent globality. We automatically obtain the structure of 
^-decomposition not only for some predefined small region, but for parameter 
space as a whole. 

2 Estimates for the number of root invariant re- 
gions in the parameter space 

Let us consider a polynomial P(s, k) with complex(for example, real) coeffi- 
cients, polynomially depending on real parameters k = (fci, . . . ,ki) G K. Com- 
plex parameters could be considered as a + jb, where a, b are real. 

This polynomial is usually considered as a characteristic polynomial of a lin- 
ear system. Vector k understood as indeterminacy or parameters of a regulator. 

We need to find such regions U C K that for all k 6 U P(s,k) will be 
continuous-time(resp. discrete-time) stable. The case of discrete-time stabil- 
ity is equivalent to continous-time one up to linear fractional transformation. 
Therefore we can consider only the continous-time case unlessotherwise is spec- 
ified. 

The /^-decomposition approach is based on study of decomposition of pa- 
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rameter space into connected regions with equal numbers of stable and unstable 
roots (root-invariant regions or regions of ^-decomposition) . 

It is clear that the border of these regions is a semialgebraic hypersurface in 
K parametrized by equalities $t(P(ju), k)) — 0, 3(P(jw, k)) — and, possibly, a 
hypersurface with an equation \a<x e& p{k)\ — 0, where (XdegP - is a leading term 
of P. 5ft and 3 here denote real and imaginary part of an expression. 

In the following theorems an upper bounds for a number of root invariant 
regions are obtained. 

Theorem 1. Let P(s, kx, . . . , k n ) be a complex(in particular, may be real) poly- 
nomial of degree t on s and of degree d on all ki together. Suppose that Jft(P(jw, kx, . . . , k n )) 
and 3(P(jw, kx, . . . , k n )), do not have any common divisors non-trivially de- 
pending on Li. 

Then the number of regions of D -decomposition is no greater than 

6(4td+M) n , 

Proof. The border of ^-decomposition in the continuous-time case is a semial- 
gebraic set in the parameter space p\ , . . . , p n . It is defined parametrically by 
the pair of equations $t(P(jw, ki, . . . , k n )) = 0, 3(P(jcj, kx, . . . , k n )) — 0, and 
an equation a t (kx, . . . , k n ) — 0. a t is a coefficient of s t in P(s, fci, . . . , fcn)f |24] 
§4.1.2). 

This semialgebraic set could be continued to a minimal algebraic variety 
containing it. Because 3l(P(joj, kx, . . . , k n ), ^s(P(ju;, kx, . . . , k n )) are coprime 
over C[w] its dimension will stay the same, therefore the number of connected 
components of complement will not be lesser. 

Consider a minimal algebraic variety X containing the set parametrized 
by $t(P(jw, kx, ■ . ■ , k n )) = 0, 3(P(jw, kx, . . . , k n )) = 0. Let g be a Groebner 
basis of an ideal generated by parametrization polynomials relative to some w- 
climinating order. Then by an elimination theorem f|10|. Ch.3 §1, Theorem 2; 
|10) . Ch.3 §1, Ex.5) X is defined by elements of g that do not contain ui. 

By Proposition 1 in ID] , Ch.3 §6 the resultant r of 5ft(P(jw, kx, ■ . ■ , k n )), 
^s(P(juj, kx, . . . , k n )) in variable u defines an algebraic set containing the border 
of ^-decomposition. 

This set could be bigger than the set defined by a basis. But because of 
coprimeness of 5R(P(joj, kx,-.-, k n )), ^s(P(jui, kx, ■ ■ ■ , k n )) up to to resultant r is 
a non-zero polynomial. Therefore the set of zeros of r divides IK in a number 
of parts that is no greater than the number of parts generated by zeros of the 
Groebner basis. 

In order to find that number we need to estimate the degree of h — rat (or 
h = r(5ft 2 (at) + 5 2 (a t ))) It can be easily shown that it is no greater that 2td+2d. 
Finally, we need to bound the number of regions. Using Theorem 2 |31] we get 
6(2deg/i)™. 

The bound for continuous-time is the same as for discrete-time one because 
linear fractional transformation does not change the degree of polynomials. □ 

It must be noted that the main result of [4] could signify that the multiplier 
6 in formulation of the theorem above could be replaced by 2. 

In 2-dimensional case we can obtain better bound, even sharp in some sense. 
To prove an upper bound for the number of root invariant regions we have to 
prove the following lemma: 
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Lemma 1. Let X be the plane real (possibly singular and reducible) affine 
algebraic curve of degree n. 

Then its complement in R 2 consists of no more than - + 2 " +2 connected com- 
ponents. 

This bound is sharp and it is reached on union of lines in general position. 

Proof. If X is a real irreducible nonsingular projective curve of degree n then, 
by Harnack inequality, ( 15 §3, [5] Ch.4 Th. 11) it divides plane to no more 
than — 2 - ) + 2 components and 2 ^ + 1 of them correspond to the 

ovals of X. 

If X is singular, then by First Harnack Theorem (}9 a , Ch.4, Th.12) real 
irreducible projective curve has no more than ( "- 1} 2 ( "- 2) -Y,h(h - 1) + 1 
branches. Here fcj are the multiplicities of singular points. 

Note that every self-intersection or self-tangency of I fragments of ovals of 
the curve in a singular point gives birth to no more than 21 new connected 
regions in a small neiborghood of a point. All other types of singularities do 
not lead to an increase of number of regions. 

Therefore if the mutiplicity of a singular point has increased by I, then the 
maximal possible number of ovals is reduced by (fcj + l)(ki + I — 1). Hence the 
number of connected components of a complement to X increases by no more 
than 21 - (k t + l)(h + 1-1), but if k it I > 1, then 21 - (k t + l)(h + 1 - 1) < 0. 
Thus singular irreducible projective curve divides the projective plane into no 
more than — — 2 ^ + 2 components. 

Let X be a reducible curve. Denote by the symbols mhn such a fact that 
m = {mi, . . . , m| m |) is a partition of n to the degrees of irreducible components 
of X. 

Let k = |{«|m* = 1}|, I = \{i\rrii > 1}\. 

We can estimate the number of regions Reg as 

Reg < max mhn (2 miem (m i -i)(m i -2)+2 + j^.^. TOiTOj . + n - k + 1) = 
= + max mh „(2|m| - k) = + max mh „(2Z + k) < 

< + max 2/+fe <„, <;,o<fc(2Z + k) = 

This estimate is obviously sharp. It is sufficient to take a reducible curve 
that consists from n lines in general position(e.g. without parallels or more than 
double intersections). 

The only thing we have to prove is the following inequality: 

(mi — l)(mi — 2) + 2 •r-^ 
Reg < max( ^ \~ m i m j + n — k + 1) 

rrti £m ~ i^j 

The first term in the left part of inequality is trivially derived from previous 
considerations. The term n — k corresponds to the transition from the projec- 
tive case to affine one. Actually, there are no more than n intersections between 
X and the infinite line, k of them are intersections between linear components 
and the infinite line, but this intersections do not add any new connected com- 
ponents. 

We can deduce from Bezout theorem that there are no more than mimj 
intersections between irreducible components of degrees m^, mj. Therefore we 
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Figure 1: Local structure self- intersection of a reducible and/or singular curve 
in a neiboghrood of self-interection point up to isotopy 

need to prove that every such intersection creates no more connected compo- 
nents than its multiplicity. 

Let us consider sufficiently small neighborhoods of points of two intersect- 
ing and non-intersecting irreducible components (see Figure 1). Those irre- 
ducible components are locally isotopic to an intersection of lines(one line if 
non-singular), therefore we need to consider only the case of lines , which is 
shown at the figure 1. 

Now we need globalize the local structure. Let us consider the following 
numeration of regions. 

Suppose that we have p irreducible components. Let us denote by N(i) a 
number of regions that the irreducible component i adds. 

Let us divide the naturals into p + 1 parts. - from 2 to 2 + N(l) — 1, from 
2 + N(l) to 2 + N(l) + N(2) - 1, from N = 2 + ££ =1 N(i) to +oo. We can 
break the plane into the finite union of neighborhoods such that each of them 
contains no more than one point of intersection of components. 

We can enumerate connected regions appearing in those neighborhoods. 

Let us take an arbitrary region that is not internal for any oval of an irre- 
ducible component. It will be the first region. 

Then we take an arbitrary neighborhood. If there are a point of intersection 
of irreducible components in it then we bypass that point counterclockwise and 
enumerate all appearing non-numbered regions using the following rule: 

1. If the component i that borders current region from counterclockwise- 
side have appeared before and all numbers between 2 + Y^w=o N(w) 2 + 
12w=o N(w) — 1 already correspond to some region, then let us numerate 
it by the first free number, that is greater than N. 

2. If, otherwise, there are free number from 2+J]^ N(w) 2+Y^ w=0 N(w) — 
1 then let us numerate our region by first free number from this range. 

Then we use the same rule to enumerate all non-enumerated regions around 
points of self-intersections of irreducible components. 
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Finally, all remaining regions will be enumerated upside-down and if this 
will not be possible - from left to right. 

We have e enumerated regions of the plane. 

Let us call a (possibly non-compact) region on the plane with border contain- 
ing more than one irreducible component a curvilinear polygon. Points (possibly 
infinite or even formal- i.e. directions in which border component changes) of 
intersections between two distinct irreducible components will be called vertices 
of a polygon. Parts of irreducible components lying betwen vertices will be 
called edges. 

Let Ri . . . R z be a curvilinear polygon. Every irreducible component form 
no more than half of all edges of R\ . . . R z . 

We have established a correspondence between the set of regions and segment 
of naturals [1, . . . , e+1). Consider it as a part of real line. Break every interval of 
form [i, i + 1) into parts corresponding to an irreducible components that appear 
in its border. The length of every subinterval in this partition is proportional 
to the number of edges formed by every component. Therefore there will be no 
element of the partition of length greater than \ — £ and lesser than 1 . 

Consider a neighborhood of any point of intersection between components. 
Let us bypass that point in any direction. During the bypass every component 
would appear to be a border for no more than 2j regions, where j is the inter- 
section multiplicity of that component in the point. From the other side the 
contribution of this component to the partition of interval [1, . . . , e + 1) in the 
point is no greater than 2j^ = j. This proves the lemma. □ 

Theorem 2. Let P(s, r,p) be a polynomial of degree t with complex(in particu- 
lar, possibly, real) coefficients. polynomially depending from two real parameters 
r, p or one complex parameter of a form r + jp. Denote by d its maximal degree 
as polynomial on variables r and p together. 

Suppose that %l(P(juj,r, p)) and Q(P(ju;, r,p)), as complex polynomials does 
not have any common divisors that nontrivially depends on uj. 

Then the number of regions of D -decomposition is no greater than 

q 2 + q + 2 
2 

with q = ltd + 2d. 

Proof. Follows from Lemma Q] and the proof of previous theorem. □ 

This bound could not be generalised onto n-dimensional case, because of lack 
of any sharp Harnack-like bounds for higher-dimensional real algebraic varieties. 
Moreover, we can easily show that the case of hyperplanes in general position 
even in 3-dimensional case do not give a maximal number of regions. 

Namely, Bihan[5 has shown that the maximal number of connected compo- 
nents of real projective algebraic surface of degree q is equivalent as q — »■ oo to 
dq 3 , where d 6 [M, ^]. Moreover [29], from Comessatti-Petrovsky-Oleinik and 
Smith-Thorn inequalities could be obtained an upper bound j^q — %q + y§9- 

But q planes could divide the space only in q +jj g+6 parts. 

Now we can write bounds for the case of polynomial matrix family. 
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Corollary 1. Let A(r,p) be a t x t-matrix, which entries are complex(possibly 
real) polynomials on real parameters r,p(may be on one complex parameter r + 
jp ) of degree no greater than d 

Suppose that %t(XA(jU)T>p)) ^(^(jw, r,p)), as complex polynomials does 
not have any common divisors that nontrivially depends on uj. 

Then the number of regions of D -decomposition is no greater than: 

q 2 + q + 2 
2 ' 

where q — 2t 2 d in continuous case and q = 2t 2 d + 2td in discrete one. 

Corollary 2. LetA(k\, . . . , k n ) be atxt-matrix, which entries are complex (possibly 
real) polynomials on real parameters k± . . .k n of degree no greater than d. 

Suppose that §l(xA(jw, k\,..., k n )) ^s{xa{ju, fci, . . . , k n )), as complex poly- 
nomials does not have any common divisors that nontrivially depends on lu. 

Then the number of regions of D- decomposition in continuous case is no 
greater than: 6(4i 2 c?)™, and no greater than 6(4£ 2 d + Atd) n in discrete one. 

3 /^-decompositions for 2-parametric families: com- 
putational study 

In order to study 2-dimensional ^-decompositions with Maple 14 we have used 
the following scheme. 

1. Convertion of family of matrices to the family of polynomials P(s,r,p), if 
needed. 

2. If we think about our family as a discrete-time object then we convert it 
to continuous-time using transformation P(s) \-t (s — l) degP P( 2 ^j). 

3. Computing an equation defining most of irreducible components of an 
algebraic curve, containing border of ^-decomposition using command 
Eliminationldeal for an ideal generated by $l(P(juj, r,p)) and r, p)) 
or eliminating u> using resultants. 

4. if P(s, r, p) has a real leading term then we add an irreducible component 
generated by its leading coefficient adeg p(f,p) else we add K 2 (adeg p(f,p))+ 
S 2 (adogp(7',p))- We get an algebraic curve containing the border of D- 
decomposition. 

5. Using PartialCylindricalAlgebraicDecomposition, we obtain a point cloud 
that contains at least one point from every region of D-decomposition. 

6. We count number of stable roots for every point from a cloud and generate 
the list of stable points. 

Let us proceed to examples. 

Example 1. Consider a discrete-time object s 6 + (r+jp)s 5 + §. Converting it to 
continuous time we obtain: p(s,r,p) = (s + l) 6 + (r+jp)(s+l) 5 (s — l) + |(s — l) 6 . 

D- decomposition of family z n + (r + jp)z n ~ 1 + a, with a > 1 has (n — 1) 2 + 1 
regions ' t T3j. In UHf . there is also given an estimate for number of root invariant 
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Figure 2: ^-decomposition and a point cloud from connected root invariant 
regions s 6 + (r + jp)s 5 + |. 

regions for families of a form a(s) + (r+jp)b(s), namely (n — 1) 2 + 2. Polynomial 
parametrization of border of D- decomposition is: 

-r - Aps 5 + 4ps - ^s 2 + 5rs 4 + 5rs 2 - rs 6 + | + ^ s 4 - §s 6 = 
-Ars - ps 6 - 3s 5 + dps 2 - p - 3s + 5ps 4 + 4rs 5 + 10s 3 = 

Explicit equation of border of D- decomposition has degree 10 : 

9216p 10 + 46080p 8 r 2 + 92160p 6 r 4 + 92160p 4 r 6 + 46080p 2 r 8 + 9216r 10 - 
-94464p 8 - 377856p 6 r 2 - 566784p 4 r 4 - 377856p 2 r 6 - 94464r 8 + 301440p 6 + 
+683136p 4 r 2 + 1051776p 2 r 4 + 276864r 6 - 309600p 4 - 619200p 2 r 2 - 309600r 4 + 
+122500p 2 + 122500r 2 - 15625 = 

By the Lemma{^ if degree of a border is 10, D -decomposition contains no more 
than 56 regions. There are 26 regions in the example considered. 

There are 56 regions if and only if the border of D -decomposition is an 
arrangement of 10 lines without parallels and triple intersections. 

At the figure 2 the D- decomposition and the point cloud are shown. There 
are no stability domains. 

Example 2. The other example flSf is a discrete-time object z 6 + (r+jp)z 5 + -^. 

It is an example of one possible fenomenon appearing in a given way of 
studying D- decomposition. Namely, an equation ^R. 2 (ae(r,p)) + $s 2 (ae(r,p)) = 
(r+ 1^) 2 = has nonzero number of roots but it does not separate root invariant 
regions - there are the same number of stable ad unstable roots from every side 
of the line defined by the equation. 
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Figure 3: The border of -D-decompotion with extraneous line that corresponds 
to the case of vanishing of the leading term and point cloud for the family 
z 6 + (r + jp)z 5 + ±. 

The border of D-decomposition(without extraneous components) is given by 
an euation of degree 10 : 

9, 216 • 10 13 p 10 + 4, 608 • lO 14 / r 2 + 9, 216 • lO 14 /?- 4 + 9, 216 • 10 15 pV+ 

+4, 608 • W u p 2 r s + 9, 216 • 10 13 r 10 + 8, 1792 • 10 13 p 8 + 3, 27168 • 10 14 p 6 r 2 + 

+4, 90752 • 10 4 Vr 4 + 3, 27168 • lO 14 ^ 2 ^ + 8, 1792 • 10 13 r 8 + 

+1, 30276416 • 10 1 V - 1, 821010752 • 10 16 A 2 + 1, 865389248 • 10 1 5 p 2 r i - 

-1, 15483584 • 10 15 r 6 + 6, 82460784 • 10 1 V + 1, 364921568 • 10 14 p 2 r 2 + 

+6, 82460784 • 10 13 r 4 + 4, 160322828658 • 10 1 5 p 2 + 4, 160322828658 ■ 10 15 r 2 - 

-3, 573226485213841 • 10 15 = 0. 

At the figure 3 the border of D-decomposition(including extraneous line) is 
shown. 

Example 3. The other example shows us how a minimal algebraic curve con- 
taining the border of D- decomposition may differ from the border of D- decomposition 
itself. 



Let us consider an output feedback problem of the form 
discrete-time system defined by matrices: 




( 0,4753 0,7579 7,9939' 
A = -0,0415 0,8905 0,7579 
\-0,0758 -0,0415 0,4753 



: ) 
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Figure 4: Minimal algebraic curve containing border of ^-decomposition and 
point cloud for an example 3. 



A minimal algebraic curve containing the border of D- decomposition con- 
sists from three intersecting ovals. But the border of D- decomposition is not an 
algebraic curve, because it does not contain part of the medium oval lying in the 
smallest one. 



Example 4. The next example is an output feedback of a form K = 
for a co ntinuous-time systems defined by matrices 
I 79 20 -30 -20\ 



A = 



-41 
167 
\33,5 



-12 
40 
9 



17 

-60 
-14,5 



13 
-38 



B 



( 0,219 0,9346^ 

0,047 0,3835 

0,6789 0,5194 

\0,6793 0,831 / 



0,0346 0,5297 0,0077 0,0668 
0,0535 0,6711 0,3834 0,4175 



This system has been considered in '281, Example 2 and '131, Example 15. 
We have found a new connected region of D -decomposition with one stable root. 
This region isaway from all other components, its border corresponds to uj > 20. 

This fact illustrates globality of our method. We don't need to prescribe any 
parameter ranges, but we do automatically get an answer for a whole space. 

Example 5. Here we give an example of polynomial family with complex coef- 
ficient for which bound from Lemma]l\ is sharp. 
Namely, let 

P(s,r,p) = (4-22j)s 2 + ((30 - 20j)r - (145 + 8j)p + 68)s+ 

+((21 - 8j)r 2 - (34 + I3j)p 2 + (40 + 54j)rp + (2 - 26j)r + (-34 + 19j> - (8 + 6j)) 



11 




Figure 5: Fragment of a minimal algebraic curve containing the border of D- 
decomposition for example 3 (right side) and fragment of Z?-dccomposition it- 
self (left side). 




Figure 6: Minimal algebraic curve containing the border of ^-decomposition 
and point cloud representing all root invariant regions. Newly discovered upper 
component separates regions with 3 and 1 stable roots. 
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Figure 7: Fragment of a minimal algebraic curve containing the border of D- 
decomposition for example 4. 
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Figure 8: The border of D-decomposition and point cloud for example 5. 



Then the border of D- decomposition could be parametrised as: 

4s 2 + (20r + 8p- 4)s - 21r 2 + AQrp + 2r - 21p 2 - 34p - 8 = 
-22s 2 - 145ps + 68s + 30rs - 54rp - 8r 2 + 19p - 13p 2 - 26r - 6 = 

j4n explicit equation of border of D -decomposition could be written as: 

1008000p 4 - 3303960p 3 r + 1782100p 2 r 2 + 978760pr 3 - 627300r 4 + 62160p 3 + 
+1366648p 2 r + 1219928pr 2 - 1200320r 3 - 1679328p 2 + 630352pr - 145904r 2 - 
-309120J3 + 310272r + 64512 

The decomposition of plane into 11 regions is shown on figure. 
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4 Conclusions 



We obtain a new technique of studying /^-decomposition using real algebraic ge- 
ometry. Method of obtaining an explicit equation for a border of ^-decomposition 
is given as well as method of constructing point clouds from all root invariant 
regions. 

An estimates for a number of root invariant regions are given in n-dimensional 
and in 2-dimensional cases. The last bound is close to the sharp one. 

This technique could be applied in construction of low order regulators as 
well as in robust stability analysis. 
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